Fluctuation theorem applied to the Nose-Hoover thermostated Lorentz gas 
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We present numerical evidence supporting the validity of the Gallavotti- Cohen Fluctuation The- 
orem applied to the driven Lorentz gas with Nose-Hoover thermostating. It is moreover argued that 
the asymptotic form of the fluctuation formula is independent of the amplitude of the driving force, 
in the limit where it is small. 
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Over the last decade, different versions of fluctuation 
formulas have been the focus of a number of publications 
in the field of non-equilibrium statistical physics. In par- 
ticular, dissipative deterministic dynamical systems with 
time-reversal symmetry have attracted some attention as 
potential candidates to model externally driven systems 
with a thermostating mechanism 0, 0, HJ . Two distinct 
results have been proposed, which, in the context of iso- 
kinetic thermostats, both characterize the fluctuations 
of entropy production. One, due to Evans and Searles 
0, is usually referred to as transient fluctuation theo- 
rem, while the other, due to Gallavotti and Cohen 
is simply known as the fluctuation theorem. The former 
addresses the fluctuations of the work done by the exter- 
nal forcing on the system, and the latter the fluctuations 
of the phase space contraction rate of non-equilibrium 
stationary states. It has been rigorously proved in the 
context of Anosov systems Q. 

To be definite, consider the externally driven periodic 
Lorentz gas with Gaussian thermostating Q. The trajec- 
tory of a particle in between elastic collisions is described 
by the equation p = E — ap, where a = E • p/p 2 is a 
reversible damping mechanism that acts so as to keep 
the kinetic energy constant, a, the phase space contrac- 
tion rate for this system is, as seen from its expression, 
equal to the work done on the particle divided by the con- 
stant temperature. Thus the work done on the particle 
is exactly compensated by the heat dissipation. In other 
words, work and heat dissipation statistics are identical 
for this system. 

Dolowschiak and Kovacs [j| recently made the obser- 
vation that work and phase space contraction rate fluc- 
tuations behave very differently for the externally driven 
Lorentz gas with Nose-Hoover thermostating. On the 
one hand, the work fluctuations, whether large or small, 
obey the Evans-Searles formula, in agreement with sim- 
ilar observations made for other systems @, Q ■ On the 
other hand, the authors observed that the phase space 
contraction rate fluctuations rapidly saturate. Moreover 
no observation of a linear regime of fluctuations in a lim- 
ited range was reported. 



The Nose-Hoover thermostated Lorentz gas on a pe- 
riodic lattice has phase space coordinates T = (q, p, £). 
Here q and p respectively denote the position and mo- 
mentum of the particle, and £ is the variable associated 
to the thermal reservoir. Between two elastic collisions, 
the dynamics is specified by the equations 



q = p, 
p = e-Cp, 

C = r rcs > 2 /(2T)-l]. 



(1) 



Here E denotes the external field, T resp the relaxation 
time of the thermostat and T the temperature |13| . 

An essential difference between the Gaussian and Nose- 
Hoover thermostated Lorentz gases is that the phase 
space of the former is compact. The fluctuation theo- 
rem Q can thus be applied to the Gaussian thermostated 
Lorentz gas as though it were Anosov 0] . For the Nose- 
Hoover thermostated Lorentz gas, the situation is differ- 
ent in that the phase space contraction rate fluctuations 
are unbounded. In that case, for large fluctuations, one 
expects a much larger probability of positive phase space 
contraction rate fluctuations. 

Nevertheless, an appropriate modification of the fluc- 
tuation theorem was given in 8J. Thus consider a time- 
reversible dissipative system with average phase space 
contraction rate (cr(T)) = <r+ > and assume this sys- 
tem verifies the chaotic hypothesis. In a language similar 
to that used in Q , the statement of the fluctuation theo- 
rem for the dimensionless contraction amplitude p is that 
there exists a finite number 1 < p* < oo, so that 
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where P T (p) denotes the probability of observing, over a 
time interval r, a fluctuation of the phase space contrac- 
tion rate (o~) T = pcr+. In particular, thinking about the 
external driving parameter in the Lorentz gas, the above 
result should be independent of its amplitude (|E| > 0). 

We argue the driven periodic Lorentz gas with Nose- 
Hoover thermostating considered in y| verifies the fluctu- 
ation relation Eq. J5J); although || correctly pointed out 
the phase space contraction rate fluctuations are bounded 
for large fluctuations, their measurements do not point 
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to the violation of Eq. (J2J) . Here we present evidence that 
the Lorentz gas with Nose- Hoover thermostating has fluc- 
tuations of the phase space contraction rate which are en- 
tirely consistent with the statement above, showing both 
saturation for large fluctuations and linear behavior for 
small fluctuations. 

As noted in the phase space contraction rate for 
this system has the form 



1 d 
Tdi 



H. 



(3) 



where cro(r) = E • p/T is the quantity relevant to the 
Evans-Searles fluctuation formula. Here H = p 2 /2 + 
TT r 2 esp £ 2 - Following the discussion in one can derive 
the distribution of a in terms of that of ctq. This was 
first done in in the framework of a Brownian particle 
dragged through water by a moving potential. In that 
case, Co has Gaussian fluctuations and one can derive 
the assymptotic form of the RHS of Eq. J5J : 



hm - log— r = f(p)a + , 

t-+oo T P T (—p) 



where 



(p-l) 2 /4, 



<p< 1, 

1 <J5<3, 

P > 3. 



(4) 



(5) 



For negative p, Hp) is odd, f(-p) = -f(p)- 

According to (5| , the fluctuations of the time averages 
of do for the Nose-Hoover Lorentz gas are Gaussian [TJj , 
which, given that the kinetic energy probability distri- 
bution is canonical in the presence of the Nose-Hoover 
thermostat, entitles us to use the result of [13, Eqs. Q- 
(0. We present in Figs. 11121 the results of numerical 
simulations on a hexagonal lattice with inter cell distance 
unity and disk radius 0.44 (consistent with the finite hori- 
zon condition). That is we take the fundamental lattice 
translation vectors to be (1,0) and (1/2,-^/3/2). Two 
different values of the external field are considered, both 
along the x-axis, E = (0.1,0) and E = (0.5,0). The 
numerical integration was performed using an algorithm 
similar to that used in Q- 

As seen from the figures, the numerical data is en- 
tirely consistent with Eq. (J2J. Indeed, as times become 
large, the small fluctuation amplitudes have a linear slope 
which approaches asymptotically the value given by the 
RHS of Eq. The data are compared to the predic- 
tion in Eqs. and show a rather good agreement, 
albeit f(p) = 2 seems to slightly overestimate the sat- 
uration level. As illustrated on the left-panel of Fig. 
this apparent discrepancy is only a finite-time correction, 
expected to be 0(1/T) [Tj. Likewise the right-panel of 
Fig. |3]shows that the slope in the linear regime of fluctua- 
tions approaches the predicted asymptotic value as aver- 
age lengths increase. Figs. 11121 moreover show the linear 
regime of fluctuations persists irrespective of the strength 
of the driving field. That is to say, the field strength is 



relevant to the asymptotic regime only through a + which 
is quadratic in the field strength, consistent with Eqs. 
(QJ-JSJ. Similar results were obtained for field values as 
low as E = (0.05, 0). The main difficulty being that the 
length of the time averages in Eq. J5J needs to be in- 
creased significantly in order for the slope to converge to 
the asymptotic value predicted by Eq. J2J). 




FIG. 1: Probability density function (above) and verification 
of Eq. (below) for different time averages, as indicated 
in the legend. The parameters are set to E = (0.1,0), T = 
1/2 and T roS p = 1. The average phase space contraction rate 
(2() = 3.25 x 10~ 3 was measured over a total of « 3 x 10 8 
collisions. The thick solid curve corresponds to the prediction 
in Eqs. Q-©. 



In summary, the driven periodic Lorentz gas with 
Nose-Hoover thermostating provides a simple example 
of a fully deterministic model whose work and phase 
space contraction rate fluctuations obey relations similar 
to those found in ^(| for the work and heat fluctuations 
of a Brownian particle in a potential well. Similar results 
were obtained in [12j for the non-equilibrium fluctuations 
of a RC circuit. This suggests an identification between 
phase space contraction rate and heat dissipation in this 
framework. The numerical evidence we presented is en- 
tirely consistent with the phase space contraction rate 
fluctuation theorem as stated in [j| . The sharp constrast 
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between this observation and the conclusion drawn by 
(similar claims can be found in Q) that the fluctua- 
tion theorem does not hold if the driving is small can be 
attributed to two reasons; (i) these authors did not con- 
sider properly the saturation of the fluctuation at larger 
p > p* , and (ii) they did not run the simulation long 
enough to observe the small but clear linear region at 
p < p* . The small deviations between our data and the 
form predicted by Eqs. Q-© should be attributed to 
finite-time corrections. We hope to further report on this 
issue in a future publication. 




FIG. 2: Same as Fig. EI for E = (0.5,0). The average phase 
space contraction rate (2£) = 7.62 x 10 -2 was measured over 
a total of ~ 3 x 10 s collisions. 
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FIG. 3: Height of saturation measured at p = 6 (for times 
T — 100, . . . , 750) (left), and slope measured about p = 
vs. 1/T of the curves shown in the bottom-panel of Fig. Qfor 
the external forcing E = (0.1,0). The vertical coordinates 
are renormalized to the asymptotic values predicted by Eq. 
©. Both curves can be linearly extrapolated to 1 as T — > oo. 
These data confirm that 0(1/T) corrections affect the finite 
time measurements. 



